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We use Bogoliubov theory to calculate the beyond mean field correction to the equation of state of 
a weakly interacting Bose gas in the presence of a tight 2D optical lattice. We show that the lattice 
induces a characteristic 3D to ID crossover in the behaviour of quantum fluctuations. Using the 
hydrodynamic theory of superfluids, we calculate the corresponding shift of the collective frequencies 
of a harmonically trapped gas. We find that this correction can be of the order of a few percent 
and hence easily measurable in current experiments. The behavior of the quantum depletion of the 
condensate is also discussed. 

PACS numbers: 



Recently, the group of Ketterle at MIT reported [lj 
the first measurements of the quantum depiction in a 
condensate of 23 Na atoms in tight optical lattices. The 
condensate fraction corresponds to the population of 
the interference peaks observed in the time-of-fiight im- 
ages, whereas the remaining diffusive background is in- 
terpreted as the quantum depletion. 

From the theoretical point of view, the quantum de- 
pletion of a weakly interacting Bose gas can be calcu- 
lated from Bogoliubov theory. Another, closely related, 
quantity is the beyond mean field correction to the equa- 
tion of state. In free space (no lattice) this is known 
as the Lee-Huang- Yang (LHY) correction 0]. In Ref.jiJ, 
Pitaevskii and Stringari suggested to observe the LHY 
correction through the measurements of the collective fre- 
quencies which are sensitive to changes in the equation 
of state. However, for tipical values of the atom den- 
sity and scattering length, this correction remains very 
small and difficult to measure. First experimental signa- 
tures of these corrections have been recently observed in 
a superfluid Fermi gas on the BEC side of the Feshbach 
resonance 0. The inclusion of the lattice is expected 
to enhance the role of correlations and therefore effects 
beyond mean-field might become visible in such configu- 
rations. Previous experimental [j| and theoretical 0,0 
attempts to investigate beyond mean field effects in the 
collective frequencies of trapped Bose gases have focused 
on quasi- ID systems, where the gas undergoes a crossover 
from a weakly interacting Bose Einstein condensate to a 
Tonks gas by decreasing the atomic density in the tube. 

In this Letter, we investigate a Bose-Einstein conden- 
sates in a tight 2D periodic potential forming a 2D ar- 
ray of weakly coupled tubes. For a fixed atom density, 
the gas is in an anisotropic 3D regime at small values 
of the laser intensity and it undergoes a dimensional 
crossover to a quasi- ID regime when the lattice depth 
is increased. We calculate the beyond-mean field correc- 
tion to the ground state energy due to quantum fluctua- 
tions along the crossover. We then include the harmonic 



trapping and apply the hydrodynamic theory of super- 
fluids to calculate the corresponding frequency shift of 
the lowest compressional mode along the axial direction, 
where atoms are free to oscillate subject to the harmonic 
potential and two-body interactions. We show that for 
values of physical parameters available in current exper- 
iments, this shift becomes large enough to be detected. 
Our results in the asymptotic ID regime are consistent 
with those obtained in Ref.0- 

Let us first consider a Bose-Einstein condensate in the 
presence of a 2D optical lattice 

V opt {x,y) = sE R [sm 2 (q B x) + sin 2 (q B y)}, (I) 

where s is the laser intensity, Eji = h 2 q B /2m is the recoil 
energy, with q B the Bragg quasi-momentum and m the 
atomic mass. The lattice period is fixed by d — ir/q B . 
Atoms are unconfined in the axial direction. The role of 
the additional harmonic trapping will be discussed later. 

The interparticle interaction is described by a s-wave 
contact potential with coupling constant g = 4irh 2 a/m, a 
being the 3D scattering length. We will discuss the situa- 
tion where the laser intensity is sufficiently large (s > 5) 
and the chemical potential ijl is small compared to the 
interband gap. We thus restrict ourselves to the lowest 
Bloch band, where the physics is governed by the ratio 
between the chemical potential and the bandwidth 8i, 
where t is the tunneling rate between neighboring wells. 
For /i <C 8£, the system retains an anisotropic 3D be- 
haviour, whereas for fi ss 8t, the system undergoes a 
dimensional crossover to a quasi-ID regime. Experimen- 
tally, this crossover can be realized by either increasing 
the laser intensity (t decreases) or increasing the atom 
density n (/i increases). 

In tight-binding approximation the states of the lowest 
Bloch band can be written in terms of Wannier functions 
as <j>kA x )<l>ky(v), where (j) k {u) = X^e itt:tI u>(M ~ and 
w(u) = exp[— u 2 /2<7 2 ]/7r 1 / 4 cr 1 / 2 is a variational gaussian 
ansatz. By minimizing the free energy functional with 
respect to a, one finds dja ~ tts 1 ^ 4 exp(— 1/4-^/s). The 
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system Hamiltonian takes the form 
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where e£ = k 2 /2m + 2t[2 — cos(k x d) — cos(k y d)] is the 
energy dispersion of the non interacting model, V is the 
volume, and g = C 2 g is an effective coupling constant, 
with C = d w4 ( u )du ~ d/y/2~Ka. 

In the following we assume that the number of atoms 
per tube is sufficiently large. Under this assumption, we 
can neglect the Mott insulator phase transition which 
would occur only for extremely large values of the laser 
intensity. Applying Bogoliubov theory to the Hamilto- 
nian ©, the energy spectrum E^ of the elementary ex- 
citations is given by E^ = e£(e k + %9 n ) 13- The ground 
state energy of the system can be obtained from the re- 
lation E = n 2 gV/2 — X)k^o( e k + 9 n ~ ^k)/2, where the 
second term corresponds to the beyond mean field correc- 
tion due to quantum fluctuations. This correction can be 
shown to be always negative (in free space this formula 
would contain an ultraviolet divergence to be cured by a 
proper renormalization of the coupling constant). 

Replacing the sum with integrals and performing the 
integration over the axial momentum k z , we find 
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FIG. 1: Left panel: scaling function f(x) [see Eq.JIJ] (solid 
line) and its asymptotic behavior (dashed line). Right panel: 
scaling function h(x) [see Eq.Q] (solid line) and its asymp- 
totic behavior (dashed lines). 



takes the asymptotic form 
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where the function f(x) is defined as 
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with 7(k) = 2 — cosk x — cosk y . Here 2F%[a, b, c, d] is 
the hypergeometric function and the integration over the 
transverse quasi-momenta is restricted to the first Bril- 
louin zone | fc^ | , | fcy | < tt. Equation Q has been inte- 
grated numerically and the result is shown in FigQ] (left 
panel). For vanishing x, f(x) saturates to the value 
4^2/3 ~ 1.89. In this limit, corresponding to 8f <C fi, 
we can neglect the Bloch dispersion and Eq. yields the 
ground state energy of a ID Bose gas 
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with linear density n\D — nd 2 and coupling constant 
9id = g/d 2 . Here L is the length of the tube. Result (0 
is in agreement with the exact Lieb-Lininger solution |9j 
of the ID model expanded in the weak coupling regime 
mgi D /h 2 ni D < 1. 

In the opposite 3D regime x ^> 1, the function 
(gj approaches the asimptotic law f(x) ~ lA3/y/x — 
16V2/15 ttx. This is shown in FigQ] (left panel) with the 
dashed line. Introducing the notation a = C 2 a, Eq.© 



where m* = h 2 /2td 2 is the effective mass associated to 
the band and a cr = —Q.2Ad^m/m*. The last term in 
the rhs of Eq.® corresponds to the generalized LHY 
correction in the presence of the optical lattice. We see 
that with respect to the free case, the LHY correction 
is magnified by the renormalization of both the coupling 
constant (5 > a) and the effective mass (to* > to). It 
is worth noticing that Eq.© is valid only for a -C 8i, 
a condition which, for large s, requires ultra-low atomic 
densities. 

The term a/d cr < in Eq.© is more subtle and 
amounts to a further renormalization of the scattering 
length due to the optical lattice. In particular, the cor- 
rect value for the scattering length for low-energy (E <C 
2-body collisions is given by 1 /a e f / = 1/5 — 1 /a cr 
In Eq.© a eff is replaced by the linear expansion 
a e ff ~ 5(1 + a/d C r). 

The equation of state u = dE/dN can be obtained 
by differentiating Eq.©. We find it convenient to write 
it as n = gn[l + k(n)], where the term proportional 
to k(n) accounts for the effects of quantum fluctua- 
tions. In the quasi-lD limit k{n) = —\/{ird 2 )y/ rag /n 
whereas in the opposite 3D regime we find k(n) = 
d/d cr + f3{m* /m)\Wn, with (3 = 32/3->/7r. 

Using the same formalism it is possible to calculate the 
quantum depletion given by N — Nq = X)k^o( e k + 9 n ~ 
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-Ek)/2£k- Taking the continuum limit, we find 
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Differently from the beyond mean field correction to the 
ground state energy the quantum depletion diverges 
for vanishing tunneling as h{x) ~ — ln(2.7a;)/v / 2- This 
signals that in the absence of tunneling there is no real 
Bose-Einstein condensation in agreement with the gen- 
eral theorems in ID In the opposite 3D regime 
x >■ 1 the function h(x) decays as 4/(3ny/2x) and from 
Eq.Q we find (N - N )/N = 8m*/(3m)(a 3 n/7r) 1 / 2 , 
which generalizes the standard 3D result in free space 

EH 

In the second part of the Letter we apply the hy- 
drodynamic theory of superfluids to investigate the ef- 
fects of quantum correlations on the collective frequen- 
cies of a trapped gas. Expanding the atom density as 
n(r,t) — n(r) + 5n(r,t), the hydrodynamic equations in 
the presence of the lattice can be written in the useful 
form [13 
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where we have introduced the notation V = 
(Vj_^/m/m*, V z ). For trapped configurations, the 
ground state density n(r) entering Eq.© has to be 
calculated by imposing the local equilibrium condition 
/io = fi(n(r)) + Vextft), where /io is the ground state 
value of the chemical potential, fixed by the proper nor- 
malization of n(r), and V ex t(r) = (u! 2 _r 2 L + ui z z 2 )/2m is 
the external trapping potential, here assumed of axial 
symmetry. 

By substituting fj, — gn[l + k(n)] into Eq.@ and re- 
taining only terms linear in k(n), we obtain 
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where iitj? (r) = [/io — V ex t(r)]/g is the Thomas-Fermi 
density profile. Equation (fTTTf) can be solved by treating 
its right hand side as a small perturbation, following the 
procedure of Ref . 0] . To this purpose one first solves the 
associated zero-th order hydrodynamic equation setting 
the rhs of Eq. (|10[1 equal to zero [ljj • The corresponding 
solution can then be used to calculate the frequency shift 
induced by the perturbation. One finds 



Sgj 9 I d 3 rdn(V 8n*)riT F dk/dnTF 
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where the integrals extend to the region where the 
Thomas- Fermi density is positive. We see from Eq. <|llf) 
that the shift in the frequency is not proportional to k(n) 
but rather to its derivative dk/dn, so that the density- 
independent term d/a cr does not contribute to the fre- 
quency shift Furthermore, the shift can be positive 
even if k(n) is negative, as it happens in the quasi ID 
regime. 

In order to observe effects beyond mean field, we focus 
on the lowest compressional mode along the axial direc- 
tion. For simplicity we assume an effective disc-shaped 
trap ui± I vW «w 2 . In this case the zero-th order dis- 
persion is given by ui = \J~?>uj z with density oscillations 
of the form 5n(r) ~ z 2 — where Z^ F = 2/i/moj 2 

is the square of the Thomas-Fermi radius along the axial 
direction. Calculating k(n) from Eq.Q and inserting it 
into Eq. (|llf) . after integration over coordinates we find 



(12) 



where n(0) — titf(O) is the density evaluated in the cen- 
ter of the trap and 
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Here x — 2t/gn(0) is the parameter controlling the 
dimensional crossover for trapped gases and G(y) 



3/2 fll 



f"{y) - 3/(y)/4y^. The function JTSJ satisfies 
y(0) = 1 and decreases monotonically as x increases, 
reaching the asymptotic law Y(x) = b/Q-xx for large val- 
ues of x. 

Equations l|12l) - (|13fl are l ne ^ e y results of this Let- 
ter. In the quasi-lD regime, corresponding to large val- 
ues of the laser intensity, Suj/lu = (21/512) y/mg/n(0)d 4 , 
showing that the shift increases by decreasing the cen- 
tral density. In the opposite 3D regime of small values 
of s, the frequency shift is instead given by Slo/lo — 
(35y / 7rm*/128m)^/a 3 n(0), and hence increases by in- 
creasing the central density. This latter equation pro- 
vides the lattice generalization of the frequency shift due 
to the LHY correction investigated in Ref.Q- 

In Fig[21 we show the calculated frequency shift as a 
function of the laser intensity for atomic samples of 23 Na 
(3D scattering length a = 2.75 nm) in a 2D optical lattice 
with period d — 297.3 nm (generated by a dye laser) for 
different values of the central density. We see that the 
frequency shift can be of the order of a few percent for 
values of s easily accessible in experiments. 

It is interesting to compare the frequency shift with the 
quantum depletion in the presence of the trap. Starting 
from Eqs Q and (JHJ , and inserting the external potential 
via a local density approximation, we find 
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s (laser intensity) 

FIG. 2: Calculated frequency shift of the axial mode for a 
bosonic cloud of 23 Na atoms in a disc trap as a function of 
the laser intensity for different values of the central density 
n(0) = 10 12 cm" 3 (solid line), 10 13 cm~ 3 (dashed line) and 
10 14 cm -3 (dashed-dotted line). The circles correspond to x = 
0.25, when the chemical potential in the center of the trap is 
equal to the bandwith [for n(0) = 10 12 cm -3 the ID regime 
is achieved for larger values of s]. Parameters used: a = 2.75 
nm, d = 297.3 nm. 



where the function Q(x) is defined as 

Q(x) = j 1 drr 2 (l - r^h (j^j ■ (15) 

For vanishing x, Q(x) diverges logarithmically as Q{x) = 
—0.26 — 7rlna;/16v2 whereas in the 3D regime x ^> 1 
we find Q(x) = \/2A\2x, yielding the known result 
(N - N )/N = (50Fm*/8m)^a 3 n(O) [l^- In Fig0we 
plot the quantum depletion evaluated from Ea. (|14|l as a 
function of the laser intensity for the same parameters 
used in Fig[3 The comparison with Fig|2] shows that the 
effects of the lattice on the quantum depletion are larger 
than on the frequency shifts. However, one should re- 
mind that the measurement of the collective frequencies 
can be obtained wich much higher precision. 

Let us finally discuss the conditions of applicability of 
our results. First, the effects of the trap are taken into 
account via local density approximation. This requires 
that the trapping frequencies to± , u> z are small compared 
to the chemical potential \x and the bandwidth 8t. In 
particular the condition lo z <C 8t ensures that the gas os- 
cillating along the tubes retains the 3D coherence. Sec- 
ond, the mean-field value ui — \/3u> z for the frequency 
of the collective oscillation is only valid for a strongly 
anisotropic trap uj±/^/m* <C u> z . For a finite value of the 
anisotropy, one should start from the more general for- 
mula obtained in Ref. valid for an arbitrary anisotropy 
in the absence of the lattice. We see from Eq.10 that 




s (laser intensity) 

FIG. 3: Calculated quantum depletion for a bosonic cloud of 
23 Na atoms in a disc trap as a function of the laser intensity 
for different values of the central density n(0) = 10 12 cm -3 
(solid line), 10 13 cm~ 3 (dashed line) and 10 14 cm~ 3 (dashed- 
dotted line). Parameters as in FigEl 



in the hydrodynamic theory the lattice enters through 
the renormalization of the effective mass m — > m* along 
the confined directions. From Ref.Q we then obtain 
(&j/w) ani = (m/Qm*)u)\/u)'l, showing that the correc- 
tion decreases by increasing the laser intensity, being pro- 
portional to m/m* . 

In conclusions, our results show that the measurement 
of the collective frequencies in the presence of 2D optical 
lattices can provide an efficient tool to investigate beyond 
mean field effects along the dimensional crossover. 
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